Anomalous Ferromagnetism of Monatomic Co Wire at the Pt(lll) Surface Step Edge 
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A first-principles investigation of the anomalous ferromagnetism of a quasi-one-dimensional Co 
chain at the Pt(lll) step edge is reported. Our calculations show that the symmetry breaking 
at the step leads to an easy magnetization axis at an odd angle of ~20° towards the Pt step, in 
agreement with experiment [P. Gambardella et at, Nature 416, 301 (2002)]. Also, the Co spin and 
orbital moments become noncollinear, even in the case of a collinear ferromagnetic spin arrangement. 
A significant enhancement of the Co orbital magnetic moment is achieved when modest electron 
correlations are treated within LSDA+[7 calculations. 

PACS numbers: 75.30.Gw, 75.75.+a, 75.10.Lp 
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Exploring magnetism in the one-dimensional (ID) 
limit has been a great challenge for many years. Only 
recently, Gambardella et al. |lj succeeded to observe 
ferromagnetism of monatomic Co wires decorating the 
Pt(997) surface step edge. By exploiting the element- 
selectivity of the x-ray magnetic circular dichroism 
(XMCD), the existence of long-range ferromagnetic order 
on Co was demonstrated below 15 K 0,0- Although the- 
oretically the Mermin- Wagner theorem |jj forbids long- 
range ID ferromagnetic order at non-zero temperatures, 
ferromagnetism in ID can be stabilized by a large mag- 
netic anisotropy energy, which creates barriers effectively 
blocking thermal fluctuations. The significance of such 
blocking mechanism was recognized earlier for the occur- 
rence of long-range magnetic order in 2D systems 0, 0| . 

The experiments of Gambardella et al. revealed novel 
magnetic properties of monatomic Co wires at Pt step 
edges. An unexpected magnetocrystalline anisotropy 
was observed: the easy magnetization axis was directed 
along a peculiar angle of +43° towards the Pt step edge 
and normal to the Co chain. The magnetocrystalline 
anisotropy energy (MAE) was estimated to be substan- 
tial, of the order of 2 meV/Co atom lj. In addition, 
a considerable enhancement of the Co orbital magnetic 
moment Mr, ~ 0.7 hb — as compared to the bulk Co Ml 
value of 0.14 hb — was deduced from XMCD experiments. 

In this paper we report a first-principles investigation 
of the anomalous ferromagnetism of a monatomic Co wire 
at the Pt(lll) surface step edge, using state-of-the-art 
electronic structure calculations. We focus on the in- 
triguing features of the quasi-lD Co wire, i.e., the easy 
axis rotated away from the (111) surface normal, the en- 
hanced Co orbital moment and huge estimated MAE. 
The key outcomes of our study are (i) the ab initio cal- 
culation of an easy axis at an odd angle rotated towards 
the Pt step edge and (ii) the prediction of an intrinsic 
noncollinearity between spin and orbital magnetic mo- 
ments of both the ferromagnetic Co wire and Pt sub- 
strate. The origin of this novel magnetic behavior, which 
is to our knowledge not present in known 2D and 3D 



itinerant ferromagnets is explained to be a consequence 
of the magnetic symmetry lowering at the surface step 
edge 0. Our calculations furthermore yield a MAE of 
the order of 4 meV/Co atom, and — using the LSDA+J7 
approach — a Co orbital moment Ml = 0.45 [Ib- 

Previously, several computational studies of the mag- 
netic properties of adatoms, clusters, and monatomic 
chains on surfaces were reported (see, e.g., 0, El E 
llljl). The calculations predict in general an enhanced 
MAE closely related to the reduced dimensionality and 
enhancement of the orbital moment. However, in all of 
these studies only transition-metal wires or adatoms on 
flat surfaces are investigated, i.e., geometries that are es- 
sentially different from the metal wire at a step edge. For 
wires, adatoms or clusters || at fiat surfaces the easy 
axis is either normal to the surface or in-line for some 
wires . 1 1 J . So far only one ab initio study of Co at a 
Pt step edge was reported, in which the XMCD spec- 
trum was computed |l2| . but the magnetic anisotropy 
was not considered. Our study focuses on the unprece- 
dented magnetic anisotropy properties observed for the 
quasi-lD Co chain. 

Methodology. We performed supercell calculations to 
model the Co chain at the Pt(lll) surface step edge. 
Supercells of various sizes were investigated. We shall 
discuss here particularly two supercells: one of small size, 
a toy model, which we name model I, and a large, realistic 
supercell, model II (see Fig. 1). Model I consists of one 
subsurface Pt layer built of 4 rows of Pt and one surface 
layer containing one row of Co atoms, two rows of Pt, 
as well as one empty row to model the step edge. Model 
II consists of a sub-subsurface and a subsurface Pt layer 
built of 6 rows of Pt atoms, while the surface step is 
modeled by 3 rows of Pt, one Co row, and two rows of 
empty Pt sites. In both supercells the vacuum is modeled 
by the equivalent of two empty Pt layers. All interatomic 
distances are adopted to be those of pure Pt. We note 
that while model II approaches the maximally treatable 
supercell size for full-potential, relativistic calculations of 
the MAE, the proportions of the experimental Co chain 
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FIG. 1: Top: schematic crystal structure of model II, used to 
represent the Co chain at the Pt(lll) surface step edge [T^ . 
Middle: Definition of the angles 8, <j>, and coordinate axes. 
The x axis is chosen parallel to [110] (along the Co wire), the 
y axis to [121] (normal to the wire), and the surface normal 
z is chosen parallel to [111]. Bottom: profile sketch with 
calculated M|/M£ values specified. 



at the Pt(997) step edge are still larger, consisting of an 
8 Pt rows wide terrace at the Pt step |2| . 

The first-principles calculations were performed using 
the relativistic full-potential linearized-augmented-plane- 
wave (FP-LAPW) method, in which the spin-orbit cou- 
pling is included in a self-consistent second-variational 
procedure . For most of the calculations the conven- 
tional (von Barth-Hedin) local spin-density approxima- 
tion (LSDA) is adopted, which is expected to be valid 
for itinerant metallic systems. In order to capture better 
the electron correlations expected for the Co 3d electrons 
in the reduced dimension also the LSDA+C7 approach, in 
the implementation of Ref . fljl has been applied. For fur- 
ther details of the calculations, see Ref. EE 

First-principles results. We first applied the conven- 
tional LSDA approach using the FP-LAPW method. To 
start with, the spin magnetization axis was chosen to be 
fixed either along the x, y, or z axis. The essential com- 
puted spin (Ms) and orbital (Ml) moments are given in 
Table IJ Table [Q reveals that the M s and M L on Co and 
Pt are noncollinear for a spin moment fixed along the y 
or z axis, but collinear when Ms is along the x axis. Non- 
collinearity of Ms and M L has been predicted previously 
|l7| for materials exhibiting a noncollinear spin magnetic 



structure, but this is to our knowledge the first observa- 
tion of such noncollinearity for a collinear, ferromagnetic 
spin configuration. To understand the noncollinearity of 
Ms and M L it is instructive to consider the magnetic 
symmetry. The symmetry operations which preserve the 
crystal symmetry are the identity E and the mirror op- 
eration a x with respect to the yz plane (see Fig. 1). Con- 
sidering now the magnetic symmetry operations, which 
are — for a total magnetic moment in the yz plane — E and 
a x R, with R the time inversion, we observe that these 
symmetry conditions impose M x —0, but M y , M z ^0 
without any restriction. Therefore there is no particu- 
lar symmetry imposed direction in the yz plane which 
would force spin and orbital moment to be parallel. In 
other words, the magnetic symmetry in the yz plane is 
the same for all magnetization directions. Along the wire 
the situation is different: the magnetic symmetry oper- 
ations E and a x , which conserve M x , force M y =M z =0 
and consequently, we must have Ml \\ Ms for a magne- 
tization along the wire. 

From Table|l]wc further observe that both the Co spin 
and orbital moment are considerabl y en hanced with re- 
spect to the values for bulk hep Co [18j, as expected for 
a dimensionality reduction leading to a more atomic-like 
configuration. The calculated Ms and M L of Co agree 
well with those of Ref. [l2|, where however only z axis 
collinear components of Ms and Ml were considered. 
The Co Ms does not change when the supercell is en- 
larged from model I to model II, but small changes in the 
orbital moments exist. Also, there is a sizable magneti- 
zation induced on the nearest neighbor Pt atoms, which 
is decreasing rapidly for the Pt atoms farther away (see 
Fig. 1). The size of supercell model II appears thus suf- 
ficient to separate the magnetic Co wires and to ensure 
the Pt magnetization decrease away from the step edge. 

Next, we turn to the salient aspect of our investigation, 
the MAE calculations. We used the so-called "magnetic 
force theorem" to compute the MAE: starting from self- 
consistent charge and spin densities calculated for the 



TABLE I: Spin and orbital magnetic moments in /j,b, calcu- 
lated for the monatomic Co wire at the Pt(lll) step edge, 
using the supercell models I and II. 
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FIG. 2: The MAE calculated for model II as a function of the 
angles 6 and <j> of the spin moment. The solid lines are fits to 
the points and are given by MAE = 4.04 - 4.44 cos 2 (0 - 18°) 
in the yz plane and 3.08 — 3.08 cos 2 in the xz plane. 



spin moment aligned along the z axis, the Ms is rotated 
over angles 9 or </> (see Fig. 1) and a single energy band 
calculation is performed for the new orientation of Ms- 
The MAE, which is defined as a directional total-energy 
difference, is computed from the change in one-electron 
energies E due to the Mg rotation, i.e., MAE = E(6, — 
E(9 = 0, <f> = 0) 0. The calculated MAE is shown in 
Fig. 2 for supercell II. 

For a rotation of M$ over an angle 9 in the xz plane 
the MAE dependence on 9 is symmetric, reflecting the 
mirror symmetry a x , with the easy axis pointing along 
the z direction and the hard axis directed along the Co 
wire. A MAE difference between the hard and easy axes 
of ^3 meV/Co (~1 meV/Co for model I) is calculated, 
exceeding by an order of magnitude the dipolar shape 
anisotropy |20|. 

For a rotation of Ms over an angle <j> in the yz plane, 
we obtain a peculiar asymmetric dependence of the MAE 
on <j> (see Fig. 2), reflecting the absence of any particu- 
lar symmetry imposed direction in the yz plane. The 
computed easy axis is rotated away from the z axis by 
18° towards the Pt step edge, in semi-quantitative agree- 
ment with the experimentally observed anomalous 43° 
easy axis Q. The calculated direction of the hard axis 
of —72° corresponds reasonably with the experimental 
value of « —50° also. A notable difference appears for 
MAE calculations adopting model I: while the calculated 
MAE(0, 9 = 0) curve is asymmetric as well, the minimum 
occurs at w —50° (not shown), thus oriented outwards 
from the Pt step. In contrast to supercell model II the 
toy model I is not valid even qualitatively for the MAE. 
Increasing the supercell size from model I to II improves 
the MAE, therefore an even better agreement with ex- 
periment can be expected for even larger supercells. The 
MAE difference between the hard and easy axes is «4.45 



meV/Co for model II (~1.6 meV/Co for model I), which 
is of the same magnitude as the experimentally estimated 
MAE of -2 meV/Co at T = 45 K ,1]. We expect a defi- 
nitely higher experimental MAE and thus an even better 
agreement with model II for T = 0K. We note that pre- 
vious studies showed the conventional LSDA theory to be 
quite successful for describin g th e uniaxial MAE of hep 
Co and CoPt bulk alloy [UlQ. Here we demonstrate 
that the LSDA also provides a reasonable explanation of 
the MAE of the Co wire at the Pt step edge. 

Although the LSDA works well for MAE calculations, 
it does not give large enough values for the orbital mo- 
ment 0. For example, the LSDA calculated Ml of hep 
Co is with 0.08 fiB only half the experimental value of 
0.14 fiB- The situation is even worse for the Co wire. 
Comparing the LSDA calculated Ml ^0.15 fis (see Ta- 
ble I) with the experimental M L of 0.68 ± 0.05 [i B Q 
indicates that the LSDA value is too small by a fac- 
tor of 4.5! Recently, the orbital-polarization correction 
(LSDA+OP) was applied to a Co wire on Pt, leading to 
a Co Ml of 0.92 fj,g E3> but this value overshoots the 
experimental data. 

To improve the Ml one needs to account for electron- 
correlation effects beyond the conventional LSDA, which 
is currently a challenging problem of ab initio relativis- 
tic energy-band theory. To estimate this effect, we use 
here the semi-model but physically transparent LSDA+C7 
method, which was shown to correct the Co Ml of 
both hep Co and CoPt alloy with a single choice of the 
Coulomb U (= 1.7 eV) and exchange J (= 0.91 eV) pa- 
rameters [21|. Using the same U and J, we compute 
M£= 0.45 fiB ror supercell II when Ms is along the z 
axis (M£ = 0.32 [Ib for model I). These values are still 
smaller than the experiment but we expect a larger Ml 
with a further increase of the supercell, due to a related 
Co d-states localization. We could of course obtain better 
agreement with experiment by another choice of U and 
J, but we prefer to use the "universal" values found in 
Ref. EH treating thus the U of metallic Co as a transfer- 
able, atom specific quantity. The Co spin moment hardly 
changes, from 2.13 to 2.18 hb when the U is included. 
We note, that in the LSDA+t/ method the MAE can be 
computed only from total energies, as it is incompatible 
with the force theorem. It will require highly accurate 
self-consistent calculations for the Ms rotated over dif- 
ferent angle 9 and 0, something which is numerically not 
a practicable approach. 

To understand how the enlargement of the Co Ml in 
the LSDA+t/ approach comes about we consider the spin 
and orbitally resolved 3d densities of states (DOS) of 
model II, which are shown in Fig. 3. The spin-resolved 
LSDA DOS reveals a substantial narrowing of the band 
width from ~6 eV for hep Co to ~4 eV for the Co wire as 
well as a moderate increase of the spin-splitting (Fig. 3a), 
as is expected for the reduced Co coordination. When the 
U is included the 3d DOS broadens somewhat and sig- 
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FIG. 3: (a) The spin-resolved 3d DOS, computed with the 
LSDA and LSDA+(7 for model II, and compared with the 
LSDA DOS of hep bulk Co. (b) The orbitally-resolved, spin- 
down DOS calculated with the LSDA and LSDA+[/ methods. 



nificant changes in the spin-resolved DOS occur. Since 
the spin-up Co (i-band is fully occupied, only changes of 
the spin-down band are essential for the Ml enhance- 
ment. The spin-down m-resolved 3d DOS of the Co wire 
is shown in Fig. 3b. A major change in the LSDA+Z7 
DOS appears as an upward shift within the |j;m = 0} 
DOS, which, however, does not contribute to Ml- The 
major contribution to the increase of Ml originates from 
an upward shift within the \l;m = —2) and a down- 
ward shift of the | J,; m = +2) DOS. Also, downward 
shifts of the ||;m = ±1} states take place, but these 
contribute only secondarily to the Ml change. Thus, we 
conclude that the Ml enhancement with the Coulomb U 
is brought about by modifications of the in-plane spin- 
down x 2 — y 2 and xy orbital densities and much less af- 
fected by changes in the out-of-plane xz, yz, and 3z 2 — r 2 
orbital densities. The in-plane orbitals are affected most 
by the missing Pt atoms at one edge side and thus most 
liable to localize. 

In conclusion, employing first-principles calculations 
we have provided a microscopic picture of the anomalous 
magnetocrystalline anisotropy of a quasi- ID Co chain at 
the Pt(lll) step edge. The essential symmetry break- 
ing at the step edge leads to noncollinear spin and or- 
bital moments as well as to an easy magnetization axis 
oriented at a peculiar angle towards the Pt step edge. 
LSDA theory is found to provide a rather good explana- 
tion of the magnetocrystalline anisotropy, yet a consider- 



able improvement of the Co orbital moment is obtained 
with LSDA+t/ calculations. 
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